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ABSTRACT 

By  examining  the  relationship  of  the  individual  waves  in 
an  ocean  wave  record  and  the  power  density  spectrum  (PDS), 
it  is  demonstrated  that  the  PDS  conventionally  obtained  using 
spectral  analysis  can  be  reasonably  approximated,  except  at 
the  very  low  frequency  end,  by  a  pseudo-power  density  spectrum 
constructed  from  the  individual  heights  and  periods  contained 
in  the  record.   The  very  low  frequencies  seen  in  the  side 
bands  of  the  PDS  are  examined  as  products  of  amplitude  modu- 
lation by  analogy  with  simulated  wave  records  and  as  an 
indication  of  the  presence  of  wave  groups.   In  addition, 
wave  groups  are  identified  and  their  properties  determined 
in  19  selected  ocean  wave  records  by  integrating  the  wave 
record  using  the  procedure  of  Funke  and  Mansard  (19  79) . 
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I.   INTRODUCTION 

The  computer  era  and  the  Fast  Fourier  Transform  (FFT) 
analysis  that  is  associated  with  it  permits  the  manipulation 
of  wave  records  in  digital  form  to  produce  a  power  density 
spectrum  (PDS)  with  great  efficiency. 

The  PDS  gives  a  presentation  of  the  distribution  of  the 
total  wave  energy  contained  in  successive  frequency  bands  or 
sometimes  in  period  bands.   In  this  way,  the  analyst  and 
interpreter  of  wave  records  has  a  process  that  describes  the 
more  or  less  random  individual  wave  heights  and  periods 
present  in  a  wave  record  with  a  deterministic  approach,  where 
the  peak  spectral  frequency  or  period  and  its  relationship  to 
the  total  energy  content  are  the  salient  features. 

However,  the  PDS  does  not  fully  describe  all  practical 
characteristics  of  the  wave  field.   For  example,  to  synthesize 
a  realistic  wave  field  in  a  laboratory  channel,  the  matching 
of  amplitude  spectra  between  the  prototype  and  the  simulated 
wave  field  is  no  longer  sufficient  for  the  reason  that 
different  wave  fields  may  have  basically  the  same  frequency 
spectrum  (Funke  and  Mansard,  1979) . 

The  particular  sequence  of  high  and  low  waves  is  an 
additional  relevant  factor.   Therefore,  attention  to  grouping 
of  waves  and  their  successive  individual  properties  must  be 
included  as  an  additional  objective  of  the  analysis  of  wave 
records . 
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Also,  it  will  be  shown  that  the  two  evident  parameters 
which  can  be  detected  by  an  observer  when  studying  sequences 
of  individual  waves,  their  heights  and  periods  (see  defini- 
tion in  Section  III. A. 2),  can  be  further  manipulated  and 
utilized  to  generate  a  so-called  pseudo-spectrum.   This 
pseudo-spectrum  may  be  considered  another  descriptor  of 
energy  in  relation  to  the  individual  waves  that  are  seen  in 
the  field. 

The  general  objectives  of  this  thesis,  therefore,  are  to 
inquire  into  the  feasibility  of  producing  a  pseudo-frequency 
spectrum  from  the  height  and  period  measured  from  the  indi- 
vidual waves  constituting  the  wave  record,  and  to  analyze 
sequences  of  heights  and  periods. 

With  regard  to  the  relation  between  spectrum  form  and 
group  properties,  it  is  also  an  objective  of  the  thesis  to 
inquire  into  their  detection  if  possible  through  considera- 
tion of  apparent  peaked  symmetrical  side  bands  occurring  in 
the  PDS  which  may  represent  the  effect  of  low-frequency 
amplitude  modulation. 

Finally,  it  is  our  intention  in  the  course  of  the  study 
presented  by  Funke  and  Mansard  (19  79)  to  examine  wave  group 
properties,  namely  their  "degree"  of  presence  in  the  wave 
record,  their  duration  and  their  variability. 
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II.   WAVE  DATA 

A.   DATA  SOURCE 

The  body  of  this  work  was  accomplished  with  data  from  the 
Point  Santa  Cruz  waverider  buoy  obtained  during  the  period 
from  June  1978  to  February  1979.   The  location,  shown  in 
Figure  1,  is  in  a  water  depth  of  60  meters.   This  is  deep  to 
intermediate  water  for  all  ordinary  gravity  waves  (sea  and 
swell)  as  defined  by  the  relative  depth. 

The  wave  buoy  was  installed  as  part  of  the  California 
Coastal  Engineering  Data  Network  (CEDN)  which  is  sponsored 
jointly  by  the  California  Department  of  Boating  and  Waterways 
and  the  U.S.  Army  Corps  of  Engineers.   Additional  funding  for 
the  wave  buoy  system  is  supplied  by  the  University  of  Cali- 
fornia Sea  Grant  Program  with  Scripps  Institution  of  Oceano- 
graphy providing  overall  direction  for  the  actual  data 
collection  program. 

The  data  from  the  buoy  received  several  times  a  day  is 
recorded  on  a  magnetic  tape  compatible  with  the  IBM  360/6  7 
computer  system.   Each  record  consists  of  10  24  discrete  points 
of  instantaneous  surface  elevation  at  1  second  interval. 
Thus  each  wave  record  has  a  duration  of  17  minutes  and  4 
seconds . 

The  individual  wave  records  are  spectrally  analyzed  by 
CEDN  and  the  results  are  published  monthly.   Tabular  data, 
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an  example  of  which  is  shown  in  Figure  2,  include: 

Significant  wave  height 

Total  energy 

Percentage  of  total  energy  in  period  bands  of  2  seconds. 

Also  a  visual  presentation  of  the  daily  energy  spectra  is 
given  so  as  to  provide  an  indication  of  the  overall  daily 
energy  pattern  for  a  one  month  period  (Figure  3) .   This 
presentation  enables  one  to  visualize  on  what  days  of  the 
month  the  energy  content  is  high,  how  the  energy  is  parti- 
tioned between  the  shorter  and  longer  period  waves,  and  how 
broad  or  narrow  is  the  spectrum. 

B.  WAVE  RECORD  ANALYSIS 

For  this  study  wave  records  were  provided  on  a  magnetic 
tape  by  the  CEDN.   Spectral  and  other  analysis  of  the  data 
were  performed  independently  of  the  CEDN  computations.   The 
published  analyses  made  by  CEDN  were  used  only  for  purposes 
of  comparison  with  certain  of  our  own  analyses.   The  wave 
records  selected  for  analysis  in  this  thesis  were  chosen  for 
their  high  total  energy  and  for  the  dominantly  unimodal 
character  of  their  spectra. 

C.  FILE  NUMBER 

For  convenience  in  discussing  individual  wave  records 
and  the  products  of  all  analyses,  file  numbers  were  assigned 
to  all  records  in  agreement  with  those  given  by  CEDN.   It 
should  be  noted  that  the  figures  are  arranged  in  a  logical 
order  but  that  they  are  not  necessarily  discussed  in  the 
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following  sections  in  the  order  in  which  they  appear.   Because 
of  the  large  number  of  printouts  contained  in  the  19  files 
analyzed  the  decision  was  made  to  include  only  one  complete 
file  in  the  thesis  for  illustration  (File  Number  228) .   The 
wave  record  for  this  file,  which  was  considered  a  good  repre- 
sentation of  a  typical  wave  field,  is  shown  in  Figure  4. 

D.   CHARACTERISTICS  OF  WAVE  RECORDS  ANALYZED 

The  significant  wave  characteristics  of  the  19  wave 
records  analyzed  are  given  in  Table  I.   The  wave  records  are 
classified  in  the  order  of  the  wave  type  as  defined  by  their 
significant  steepness  (after  Thompson  and  Reynolds,  1978), 
where  steepness  is  defined  as 

^  HS 


g  Tp2 

by  analogy  with  linear  theory,  Hs  is  significant  wave  height 

and  Tp  is  peak  spectral  period.   The  defining  steepness 

limits  adopted  were  the  following: 

Sea  1/12  -  1/40 

Young  swell  1/40  -  1/100 

Moderate  swell  1/100  -  1/250 

Old  swell  >l/250 
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III.   INDIVIDUAL  WAVE  ANALYSIS 

This  chapter  considers  the  question  of  whether  the 
frequency  spectrum  of  a  wave  record  can  be  approximately 
constructed  from  the  H  and  T  measures  of  the  individual  waves 
constituting  the  record.   If  this  is  shown  to  be  feasible,  it 
is  felt  that  it  will  permit  a  fuller  interpretation  of  the 
frequency  spectrum  in  terms  of  the  waves  composing  it. 

The  approach  taken  in  addressing  this  question  was  to 
produce  a  true  spectrum  from  a  given  wave  record  using  the 
FFT  technique  and  also  a  so-called  pseudo- spectrum  using  the 
H  and  T  measures  from  the  same  record. 

For  various  reasons,  spectra  are  difficult  to  compare. 
Accordingly,  both  spectra  were  converted  into  cumulative 
energy  distributions  and  these  were  then  compared.   These 
procedures  and  the  principal  results  obtained  are  described 
below. 

A.   SPECTRAL  ANALYSIS 
1 .   True  Spectrum1 


1  Although  the  values  are  "spectral  estimates"  of  the 
"true  spectrum",  with  uncertainty  described  by  the  Chi-squared 
probability  density  function  with  eight  degrees  of  freedom, 
the  name  true  spectrum  was  chosen  as  contrasted  to  the  "pseudo- 
spectrum"  derived  from  the  individual  H^  and  Tj_. 
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a.  Power  Density  Spectrum  (PDS) 

With  a  basic  program  that  utilizes  the  Inter- 
national Mathematical  and  Statistical  Library  subroutine 
RHARM,  the  power  density  was  found  via  the  Fast  Fourier 
Transform  performed  on  each  record  of  length  10  24  points 
with  one  second  interval. 

The  data  were  divided  into  four  "windows"  which 
results  in  spectral  estimates  distributed  according  to  a 
Chi-squared  density  with  eight  degrees  of  freedom,  while 
maintaining  sufficient  resolution  to  detect  the  low-frequency 
components  for  periods  as  long  as  four  minutes.   Thus  the 
width  of  each  window  is  256  seconds,  the  maximum  period  that 
can  be  detected.   Consequently,  the  resolution  is  1/256  Hz 
or  0.0039  Hz. 

It  is  worthwhile  to  note  that  for  more  "windows" 
of  shorter  duration,  say  8  with  128  seconds  width,  the  number 
of  frequency  components  in  the  PDS  decreases.   This  increases 
the  stability  of  the  spectral  estimates,  since  in  the  fewer 
but  broader  spectral  bands  the  number  of  degrees  of  freedom 
is  increased.   For  comparison,  Figures  10  and  11  show  PDS 
produced  for  the  same  wave  record  (File  Number  228)  using 
four  and  eight  windows,  respectively  (note  that  the  energy 
density  scales  differ) . 

b.  True  Energy  Histogram 

Also  using  a  class  interval  of  Af  =  1/256  Hz  a 
histogram  of  the  distribution  of  energy  in  percentage  of 
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the  total  was  computed  (Figure  12)  from  the  PDS  (Figure  10) . 
This  true  energy  histogram  represents  the  distribution  of 
energy  centered  about  the  same  frequencies  that  define  the 
PDS.   In  a  sense  each  value  represents  the  area  of  a  "slice" 
resulting  from  the  sum  of  two  pieces  with  base  width  Af/2 
extending  on  each  side  of  the  central  frequency. 

For  comparison,  an  energy  histogram  with  class 
interval  1/128  Hz  produced  from  the  PDS  obtained  using  eight 
windows  (Figure  11)  is  shown  in  Figure  13.   In  both  figures, 
these  histograms  are  labelled  true  energy.   The  values  are 
connected  by  a  line  as  an  aid  to  visual  inspection, 
c.   Cumulative  True  Energy  Distribution 

Early  trials  showed  that  comparison  of  the  true 
energy  with  our  pseudo-energy,  not  described  yet,  was  better 
done  with  a  cumulative  curve,  summed  from  low  to  high 
frequency.' 

The  normalized  cumulative  energy  distribution 
computed  from  the  energy  histogram  in  Figure  12  is  shown  in 
Figure  15  (labelled  true  spectrum),  with  ordinates  of 
percentage  of  the  total  energy. 

A  cumulative  energy  description  in  energy  units 
of  length  squared  was  also  produced.   Figure  14  presents 
this  curve  also  labelled  true  spectrum. 
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2.   Pseudo -Spectrum 

a.  Determination  of  Hj_  and  Tj_ 

Using  the  zero-upcrossing  method,  the  period, 
Tj_,  and  the  height,  Hj_,  of  successive  individual  waves  were 
measured  in  a  wave  record.   The  still  water  level  was 
obtained  from  the  mean  of  the  10  24  instantaneous  surface 
elevation  digital  points  (the  zero  level  reference) .   The 
individual  period,  Ti,  was  taken  as  the  time  difference 
between  two  successive  zero-upcross  points.   The  individual 
wave  height,  Hj_,  was  taken  as  the  crest  to  trough  distance 
where  the  crest  and  trough  were  represented  by  the  maximum 
positive  and  minimum  negative  values,  respectively,  within 
each  pair  of  zero-upcross  points. 

b.  Pseudo-Energy  Histogram 

The  energy  histogram  that  results  from  a  power 
density  spectrum  (PDS)  is  a'  histogram  of  energy  versus 
frequency,  where  the  energy  is  interpreted  as  energy  per 
unit  area  of  sea  surface  per  unit  weight,  with  dimensions  of 
length  squared.   This  energy  is  associated  with  the  central 
frequency  of  the  band. 

It  is  also  possible  to  treat  a  sequence  of  waves 
so  that  individual  Hj_  and  T£  are  presented  in  terms  of  Hj_2 
and  1/Tj_,  the  same  dimensions  as  above.   It  is  logical  that 
the  average  of  the  Hj_2  in  the  wave  record  should  represent 
a  parameter  that  is  directly  related  to  the  energy  per  unit 
area  because  those  are  the  waves  in  that  space. 
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Thus  if  each  Hj_2  divided  by  the  total  number  of 
waves,  N,  in  the  record  lends  its  contribution  to  the  average 
total  energy  per  unit  area,  we  can  use  it  as  an  energy  para- 
meter for  waves  within  a  class  interval  as  well.   Then,  the 
same  class  intervals  previously  used  to  produce  the  true 
histogram  can  produce  a  pseudo-energy  histogram. 

Thus,  in  our  approach,  the  individual  wave  period 
was  converted  to  its  reciprocal.   This  parameter,  with  the 
dimensions  of  frequency,  was  then  organized  into  class  groups. 

Hi2 
Then  the  summation  of  in  each  class  interval  was  computed 

and  this  value  was  found  in  percentage  of  the  total  sum  in 
all  classes.   With  the  central  frequency  taken  as  reference, 
a  plot  of  the  summed  values  in  each  class  interval  yielded  a 
histogram.   Figures  12  and  13  present  histograms,  labelled 
pseudo-energy,  plotted  using  class  intervals  of  1/256  Hz  and 
1/128  Hz,  respectively.   These  are  superimposed  for  comparison 
on  the  true  energy  histograms  derived  from  the  PDS .   Both 
sets  of  curves  are  produced  from  the  same  wave  record  shown 
in  Figure  4  (File  Number  228) .   A  comparison  of  the  two 
curves  is  discussed  later. 

c.   Cumulative  Pseudo-Energy  Distribution 

In  order  to  compare  more  objectively  the  pseudo- 
energy  histogram  derived  from  H^  and  T^  values  with  the  true 
energy  histogram  derived  from  the  PDS,  the  histograms  were 
transformed  into  cumulative  energy  distributions.   This  was 
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done  by  cumulating  the  energy  in  each  class  interval  (frequency 
band)  from  low  to  high  frequency  in  the  same  manner  as  with 
the  true  energy  histogram.   Figure  15  provides  an  example  of 
the  cumulative  pseudo-energy  distribution  in  percentage  for 
File  Number  228,  derived  from  the  pseudo-energy  histogram 
with  a  1/256  Hz  class  interval  (Figure  12) . 

With  the  objective  of  standardizing  the  total 
energy  in  the  pseudo-spectrum  to  the  total  energy  in  the  true 
spectrum  in  energy  units  of  length  squared,  an  attempt  was 
made  to  obtain  a  conversion  parameter.   This  was  accomplished 
by  dividing  the  summation  of  Hj_2  by  the  number  of  waves  in 
the  record,  N,  and  by  the  variance,  V,  obtained  from  statis- 
tical analysis, 

N 

h  =  *i2 

N  i=1   x 


V 

The  interesting  result  is  that  a  mean  value  of  6.70  with  a 
standard  deviation  of  0.37  was  obtained  for  the  19  records 
analyzed.   These  individual  values  are  presented  in  Table  II. 
The  reason  that  the  total  energy  factor  appears  to  be  a 
constant  and  to  have  the  value  6.70  is  not  clear. 

In  view  of  this  apparent  constant  it  was  decided 


to  apply  a  proportionality  factor  of  6.7  to  convert  our 


Hi 
individual  energy  parameter,  —zz —  ,  to  a  new  parameter 


Hj2 

6.7xN 
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Following  this,  each  individual  wave  was  put  in  ascending 

order  by  its  apparent  frequency  (the  inverse  of  its  period) 

Hi2 

and  each  individual  was  summed  in  a  cumulative  way. 

6.7xN  * 

The  results  of  File  Number  228  and  of  four  typical  records 

chosen  out  of  19  records  are  presented  in  Figures  14,  32,  42, 

5  2,  and  6  2  in  comparison  with  the  cumulative  true  energy 

distribution. 

3.   Comparison  of  Pseudo  and  True  Spectra 

Since  the  presentation  of  the  cumulative  energy 
distribution  in  energy  units  depends  on  the  proportionality 
factor  chosen,  the  cumulative  curves  were  compared  in 
percentage  (see  examples  in  Figures  15,  32A,  42A,  52A  and  6  2A) 

The  significant  points  that  result  from  comparing  the 
two  cumulative  energy  curves  derived  from  the  19  wave  records 
analyzed  are  presented  below: 

The  two  curves  agree  quite  closely  about  the  central 
part  of  the  distribution,  but  differ  most  toward  both  ends. 
Thus  the  slope  of  the  central  portion  of  the  pseudo-spectrum 
identifies  the  true  spectrum  peak  and  also  characterizes  the 
proper  spectral  width. 

The  true  spectrum  in  all  of  the  records  analyzed 
identifies  a  significant  amount  of  energy  in  frequencies 
lower  than  the  frequency  equivalent  to  the  reciprocal  of  the 
maximum  wave  period  detected  by  the  zero-upcrossing  method. 
At  this  frequency  it  was  found  that  the  differences  in 
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percentage  of  cumulated  energy  between  the  true  and  pseudo- 
spectra  among  the  19  wave  records  ranges  from  1.75  to  18. 
Thus  the  lower  frequencies  of  the  true  energy  distribution 
are  poorly  represented  in  the  pseudo-energy  distribution. 

It  is  also  evident  that  the  pseudo-spectrum  differs 
substantially  from  the  true  spectrum,  in  the  sense  that  it  is 
less  smooth.   This  means  that  the  same  quantity  of  energy  is 
described  with  differently  weighted  components  by  the  pseudo- 
spectrum. 

B.   WAVE  HEIGHT,  PERIOD,  AND  STEEPNESS 

As  already  described,  the  production  of  the  pseudo- 
cumulative  spectrum  starts  with  the  computation  of  the 
successive  wave  heights  and  periods.   It  was  of  some  interest 
for  the  purpose  of  visually  comparing  wave  heights  and  periods 
with  both  the  true  and  the  pseudo- spectra  produced,  and  also 
in  regard  to  the  examination  of  wave  groups  occurring  in  the 
wave  record  (described  later) ,  to  print  out  in  graphical  form 
the  sequence  of  Hj_  and  T±    occurring  in  the  record.   From  our 
19  files  we  present  the  sequences  of  heights  and  periods  for 
five  wave  records  in  Figures  5,  6  (File  Number  228),  24,  25 
(File  Number  579),  34,  35  (File  Number  600),  44,  45  (File 
Number  446),  and  54,  55  (File  Number  677). 

In  addition,  the  individual  wave  steepness  derived  from 
Hj_  and  Tj_  was  also  computed  for  each  wave  and  the  sequence  of 
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steepness  values  was  also  printed  out  graphically  for  each 
wave  record  (examples  are  shown  in  Figures  7,  26,  36,  46,  and 
56)  . 

The  individual  wave  steepness,  defined  by  analogy  with 
the  linear  wave  theory,  is  given  by 

H.j        2ttHt 


Li       g   Ti2 

We  are  aware  that  in  real  water  waves  the  length  of  an 
individual  wave  may  be  different  from  the  length  given  by 
linear  theory.   Nevertheless,  we  reasonably  assume  that  the 
quantity  Sj_  defined  here  is  proportional  to  the  height-to- 
length  ratio  in  real  waves. 

Finally,  we  present  in  Figures  8,  27,  37,  47,  and  57 
graphs  of  H-  versus  Tj_;  these  can  be  interpreted  as  a 
representation  of  the  individual  wave  steepness  distribution 
in  the  record.   Also  plotted  on  each  graph  is  the  significant 
wave  steepness, 

27TH,, 
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S  m   2 
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which  characterizes  the  wave  record. 

1 .   Sequences  of  Individual  Height  and  Period  (Hj_,  T-J 
The  following  is  a  summary  of  the  most  significant 
points  that  were  found  from  examination  of  the  sequence 
graphs  of  Hj_  and  Tj_: 
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The  sequence  of  wave  heights  gives  a  first  indication 
about  the  presence  of  wave  groups  and  permits  one  to  count 
the  number  of  waves  in  a  group  (see  Figures  35,  45,  and  55) . 
It  is  interesting  to  note  that  the  number  of  successive  waves 
that  exceeds  the  significant  wave  height  is  normally  one  to 
two  and  rarely  three  to  four  or  more  waves.   It  would  be 
interesting  to  compare  these  observations  with  the  theory  of, 
e.g.,  Ewing  (19  73)  on  the  length  of  runs  of  high  waves  as 
a  function  of  spectral  bandwidth,  but  that  is  beyond  the 
scope  of  this  thesis. 

In  the  sequence  of  periods,  the  spectral  peak  is 
consistently  higher  than  the  mean  period,  as  expected. 
However,  in  records  that  can  be  characterized  as  "sea",  the 
difference  is  small  (Figure  24) ,  and  the  sequence  of  periods 
appears  to  be  random;  whereas,  when  looking  at  records  of 
young  and  moderate  swell,  such  as  shown  in  Figures  44  and  54, 
well-defined  wave  groups  can  actually  be  distinguished  by 
their  higher  values  and  patterned  sequence. 
2.   Sequence  of  Individual  Steepness,  S^ 

It  was  found  that  records  classified  as  "sea"  (Table 
I)  display  more  uniform  values  of  steepness  (Figure  26) , 
whereas  steepness  values  in  young  and  moderate  swell  have 
higher  variability  (see  Figures  36  and  46;  note  the  different 
vertical  scales) .   Figure  56,  as  expected  with  old  swell, 
presents  a  preponderantly  constant  steepness.   All  sequence 
graphs  of  Sj_  characteristically  show  steepness  spikes. 
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3.   Individual  Height  Versus  Period 

Each  plotted  point  in  these  graphs  represents  the 
steepness  S-  of  an  individual  wave  in  the  wave  record.   The 
points  are  seen,  in  most  of  the  graphs,  to  fall  within  a 
generally  well-defined  steepness  envelope.   The  envelope  is 
associated  with  large  steepness  values  for  sea  and  young 
swell  (Figures  27  and  37)  and  low  values  for  old  swell 
(Figure  57) ,  as  expected.   In  most  graphs  the  width  of  the 
envelope,  representing  the  scatter  of  the  Sj_  values,  is 
relatively  narrow  (Figure  27)  but  in  several  it  is  broad 
(Figure  66)  . 

In  some  graphs  the  data  points  show  a  pattern  suggesting 
more  than  one  steepness  envelope.   In  Figure  8,  for  example, 
the  data  points  associated  with  the  highest  and  steepest 
waves  are  contained  within  the  steepness  values  of  1/40  and 
1/100,  whereas  a  set  of  lower,  less  steep  waves  lies  within 
the  envelope  bounded  approximately  by  the  steepness  values 
of  1/100  and  1/250.   In  Figure  57,  the  higher  waves  have  the 
lower  steepness,  their  envelope  being  centered  about  Sj_  =  1/250, 
whereas,  a  second  set  of  waves  having  low  heights  but  high 
steepness  appear  to  occupy  an  envelope  contained  between  the 
steepness  values  of  about  1/40  to  1/50 . 

The  significant  steepness  of  the  wave  record  which 
is  also  shown  on  each  graph  is  generally  representative  of 
the  waves  of  greatest  height  and  period  in  the  record,  as 
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may  be  expected.   However,  it  was  found  that  for  a  few  records, 
significant  steepness  is  not  a  good  representation  of  the 
waves  observed  (Figure  66) .   This  situation  may  be  due  to  the 
presence  of  a  young  sea  superimposed  on  an  old  swell,  as 
Figure  6  3  suggests. 

C.  AUTOCORRELATION  OF  SURFACE  ELEVATION 

The  application  of  the  autocorrelation  function  to  a  wave 
record  and  the  analysis  of  its  periodicity  is,  as  is  known, 
another  way  of  finding  the  main  spectral  frequency  of  the 
record  (see  Figures  9,  28,  38,  48,  and  58) .   However,  this 
process  can  also  provide  an  indication  of  the  presence  of 
wave  groups. 

As  an  example,  compare  the  autocorrelation  of  a  sea  record 
(Figure  28)  with  an  old  swell  record  (Figure  58) .   The  small 
exponential  decay  in  the  swell  correlogram  is  in  correspondence 
with  the  form  of  the  envelope  of  the  wave  record  (Figure  53) , 
and  is  a  distinct  "signature"  of  the  presence  of  well-defined 
wave-groups,  in  contrast  with  the  rapid  decay  observed  in  the 
sea  type  wave  record. 

D.  COMPARISON  OF  TOTAL  ENERGY  PARAMETERS 

In  the  analysis  of  wave  records  to  obtain  the  cumulative 
true-energy  and  cumulative  pseudo-energy  distributions  several 
estimates  of  the  total  energy  were  obtained.   Later  another 
function  designated  SIWEH,  but  not  described  yet,  was 
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obtained  whose  mean  value  is  also  an  indication  of  the  total 
energy.   With  the  objective  of  comparison,  Table  III  presents 
for  the  19  wave  records  analyzed  the  total  energy  computation 
obtained  using  the  following  procedures: 

-  Statistical  variance 

-  Integration  of  the  PDS 

-  Mean  of  the  Smooth  Instantaneous  Wave  Energy  History  or 
STWEH  (Section  IV-B) 

-  Cumulative  pseudo-energy  using  a  factor  of  proportionality 
of  6.7  (Section  III-A-2(c)} 

-  Total  energy  presented  by  CEDN 

From  this  table  it  can  be  said  that  the  first  three  of 
these  methods  of  computing  the  variance  from  a  wave  record 
present  almost  equal  results,  the  differences  being  due  to 
the  differences  in  the  methods  of  numerical  computation.   It 
is  of  interest  that  the  cumulative  pseudo-energy  computation 
using  a  factor  of  6.7  gives  a  value  for  the  total  energy  with 
an  average  error  of  plus  or  minus  5.3  percent  in  relation  to 
the  variance. 
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IV.   WAVE  GROUP  ANALYSIS 

The  question  examined  here  about  'what  process  may  be 
interpreted  as  operating  in  the  peak  spectral  frequency  was 
due  to  the  graphs  showing  the  sequence  of  Hj_  and  T^  that  were 
produced  for  the  identification  of  individual  waves  and  also 
from  certain  characteristics  of  the  PDS . 

Both  sequence  graphs  show  high  variability  and  if  the 
sequence  of  wave  height  as  manifested  by  wave  groups  can  be 
interpreted  as  the  consequence  of  amplitude  modulation,  the 
sequence  of  periods  suggests  that  some  kinds  of  frequency 
modulation  may  also  be  present.   We  have  to  say  to  the  reader 
that  we  did  not  find  conclusive  evidence  for  a  simple  and 
easily  identifiable  frequency  modulation  in  our  wave  records 
examined,  but  the  question  is  raised  because  it  seems  that 
more  could  be  done  in  describing  the  wave  record  through  a 
more  deterministic  approach  than  has  been  taken  before.   Thus 
we  are  assuming  that  if  the  PDS  of  real  wave  records  show 
peaks  on  either  side  of  the  central  peak,  this  should  be  due 
to  amplitude  modulation  instead  of  frequency  modulation. 


2  The  possibility  should  not  be  ignored  that  secondary 
peaks  can  arise  as  artifacts  of  the  analysis  procedure  due 
to  "leakage" .   If  so,  they  appear  at  intervals  related  to 
the  resolution  Af. 
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It  was  felt  that  a  simple  mathematical  approach  to  modu- 
lation phenomena  will  help  in  the  interpretation  of  narrow 
band  PDS  and  in  probing  the  nature  of  wave  groups  in  real 
wave  records.   Accordingly,  this  is  presented  in  some  'detail 
in  the  next  three  sections. 

A.   MODULATION  AND  WAVE  GROUPS 
1 .   Amplitude  Modulation 

In  this  type  of  modulation  the  amplitude  of  the 
carrier  wave,  or  carrier,  is  changed  in  a  proportional  way 
to  the  instantaneous  value  of  the  modulating  signal.   We  can 
represent  this  concept  by: 
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Wave  modulated  in  amplitude 
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For  simplification  we  assume  that  the  modulating  wave 
and  the  carrier  are  sine  waves.   In  this  way  we  can  produce 
something  with  this  form: 


ec  =  Ec  sin  2TTfct       (1) 


(2) 


6  =(Ec  +  em)  sin  2Trfct  (3> 


We   can  manipulate   expression    [3]    in   this    form: 
e      =       (Ec   +    em)     sin    2i\fQt 


=      (Ec  +   Em    sin   2irf_t)    sin    2iTf    t 


•c   T  ^m 


m' 


=      Ec    sin    2TTfct   +   Em   sin    2iTfct    sin    27Tfmt 


Keeping    in  mind   that 


sin   a    sin   b      = 


cos    (a-b)  cos    (a+b) 


we   get: 


m 


e      =      Er    sin    2iTf_t   +   -=-  cos    2ir(f_   -    f_) 


m' 


jm 


-     -f    COS     2TT(fc    4-     fjt 
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Now  if  we  define  as  an  index  of  modulation  the  value  m  = 


.3 
Eo   ' 


we   get: 


mEc 
e      =      Ec    sin    277fct   +  — 5—   cos    27r(fc    -    fm)  t 

mE„ 


T"  cos  2TT(fc  +  fm}t 


This  is  the  general  expression  for  a  wave  train  modulated  in 
amplitude  by  a  single  frequency. 

What  can  be  learned  from  this  is  that  the  spectrum 
of  such  a  modulated  wave  should  present  three  distinct  peaks 
at  the  frequencies  fc,  (fc  +  fm) ,  and  ( f c  -  fm) . 

In  Figures  68  (lower  graph)  and  69,  we  have  a  simple 
wave  form  and  its  spectrum,  with  Ec  =  100  and  fc  =  0.1  Hz 
(T  =  10  sec)  modulated  by  an  Em  =  80  with  f  =  0.016  Hz  (T  = 
60  sec) .   The  spectrum  for  this  simple  case  can  be  interpreted 
in  this  way: 
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The  diagram  shows  a  central  peak,  two  side  bands  equally 
spaced  from  the  central  peak,  and  a  low  frequency  component. 

As  far  as  we  can  understand,  the  modulating  wave  in 
true  wave  records  is  not  as  simple  as  presented,  but  can  be 
visualized  by  a  complex  of  frequency  components. 

In  a  visual  sense,  it  seems  possible  that  the  charac- 
terization of  the  low-frequencies,  fm,  may  be  achieved  by 
extrapolating  to  the  following  schematic  concept: 
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In  this  representation  it  can  be  understood  that  the 
existence  of  symmetrical  side  bands  can  be  an  indication  of 
a  well-defined  low-frequency  modulating  signal,  and  thus  the 
presence  of  symmetrical  side  bands  in  the  PDS  may  be  expected 
to  indicate  the  occurrence  of  wave  groups.   So  the  wave  group 
frequency  might  be  deduced  from  the  frequency  separation 
between  the  side  bands  and  the  spectral  peak. 
2 .   Frequency  or  Phase  Modulation 

In  frequency  modulation  the  instantaneous  value  of 
the  carrier  frequency  is  changed  by  the  modulating  signal, 
but  the  amplitude  of  the  carrier  is  kept  constant.  Simple 
modulation  of  sinusoidal  waves  is  illustrated  as  follows: 


Modul ator 
signal 


Carrier 


Wave 
modul ated 

in 

frequency 
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Here  the  instantaneous  difference  between  the  carrier 
frequency  and  the  frequency  of  the  modulated  signal  is 
proportional  to  the  instantaneous  amplitude  of  the  modulator 
signal,  but  is  independent  of  the  frequency  of  the  modulator 
signal.   However,  the  number  of  times  per  second  that  the 
frequency  changes  around  the  central  frequency  is  the  frequency 
of  the  modulator  signal.   Thus  the  maximum  excursion  of  the 
frequency  from  the  central  or  peak  frequency,  Af  is: 

Af   =   K  Em 
where:   Af   =  maximum  change  in  frequency 

Em  =   amplitude  of  the  modulating  signal 

K    =   constant  of  proportionality 
To  understand  the  relationship  between  side  bands  and  carrier, 
let  us  consider  that  the  initial  carrier  can  be  written  as: 

e   =   E  sin  8(t) 
where:   e    =   instantaneous  amplitude 

E   =  maximum  amplitude 
9(t)   =   angular  value  at  time  ( t) 

The  instantaneous  angular  velocity  is  by  definition 
the  rate  of  change  of  9(t)  with  time, 

36(t)       n  c 
UjL   =   — ^ =   2irfi 

Thus   we   can    say    that   our    carrier   with    constant 
frequency   fc   has   ooc   =    2i\iQ   and   our    6(t)    =    ooct   +^   ,    where    ^ 
is    the    initial    phase   angle   at   t   =    0 . 
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Here,  as  the  frequency  is  constant 

99(t) 

wi  -  —Ft —  =  wc 

which  corresponds  to  the  fc  (i.e.,  coc  =  2irf  )  . 

Now,  according  to  the  previous  definition,  a  sinu- 
soidal signal  modulated  in  frequency  by  another  sinusoidal 
signal  has  an  instantaneous  value  of  frequency,  expressed  by: 

fj_   =   fc  +  Af  cos  (umt  , 
where:   Af   =   K  Em 

and  the  instantaneous  angular  velocity  should  be  (multiplying 
both  terms  by  2tt) 

^i   =   ^c  +  ^Tr  Af  cos  °°mt 

where:   oo^   =   instantaneous  angular  velocity 
coc   =   angular  velocity  of  the  carrier 
u)m   =   angular  velocity  of  the  modulating  signal 
Af   =   maximum  change  in  frequency 

36(t) 


Now  if      to 


r       3t 


then    6(t)       /    a>i(t)  dt 


which  gives      8  ( t)   =   to„t  +  sin  co  t 

c      um        m 

therefore  the  amplitude  of  our  initial  wave  when  under  the 
influence  of  the  modulating  signal  can  be  expressed  by: 


e   =  E  sin  (^ct  +  sin  wmt) 

wm 
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Based  on  this  expression,  the  concept  of  index  of 
modulation  in  frequency  is  defined  by 


m 


Af        maximum  change  in  frequency 
f       fm       frequency  of  modulating  signal 


Keeping  in  mind  that  Af   =   K  Em  ,  which  means  that 
Af  is  proportional  to  the  amplitude  of  the  modulating  signal, 
the  index  of  modulation,  assuming  a  constant  Af,  is  inversely 
proportional  to  the  frequency  of  the  modulating  signal. 

The  side  bands  originated  by  frequency  modulation, 
instead  of  a  simple  pair  separated  from  the  fc  by  the  value 
of  f   as  in  amplitude  modulation,  are  separated  by  Af 's  that 
are  multiples  of  the  modulating  signal.   The  number  of  extra 
side  bands  depends  on  the  relation  between  the  frequency  of 
the  modulating  signal  and  Af.   In  a  sense  this  depends  on 

m£   =   — = —   (index  of  modulation  in  frequency)  . 
rm 

It  can  therefore  be  shown  that  the  development  of  the 

wave  form  with  frequency  modulated  by  a  sinusoidal  wave,  as 

expressed,  can  be  presented  in  the  following  manner: 

e   =   E  /   (nif)  sin   id  t 
Jo 

E    +     /        (mf)        [sin    (wc    +   oo    )  t   -    sin    (ojc    -    u>    )  t] 
E    +     /        (mf)        [sin    (coc   +    2wm)  t   +    sin    Uc    -    2a>m)  t] 


+ 


+   E     /  (m.p)        [sin    (u„   +   nco    )t   +    sin    (to      -    noom)t] 

t  J  .    >  £  cm  cm 

^  (n) 
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where   /    (m^)  are  the  Bessel  functions  of  order   n  with 

argument  m^ . 

The  objective  for  this  approach  is  not  only  to  show 
that  a  typical  frequency  modulation  has  these  kinds  of 
components  from  a  spectral  point  of  view,  but  also  to  take 
note  of  the  fact  that  for  an  it\£   =   0.6  the  lateral  bands  of 
second  and  higher  order  are  approximately  zero  because  the 
Bessel  function  with  that  argument  is  almost  zero.   This  means 
that  only  the  fundamental  and  the  first  harmonic  are  important 
in  regard  to  the  spectrum  of  such  a  signal,  so  that  modulation 
of  amplitude  and  frequency  from  the  spectral  point  of  view 
look  the  same . 

3 .   Modulation  in  Real  Waves 

To  exemplify  the  special  equivalence  described  in  the 
previous  section  and  with  the  objective  to  better  interpret 
real  waves,  the  reader  is  asked  to  compare  Figures  69  and  73, 
derived  from  the  wave  records  shown  in  Figures  6  8  and  71, 
respectively  (lower  graphs)  .   As  stated,  the  PDS  in  Figure  69 
was  produced  from  amplitude  modulation  with  f    =   1/10  Hz 
and  fm   =   1/60  Hz  having  E    =   100  and  Em   =   80,  whereas 
the  PDS  in  Figure  73  was  produced  from  frequency  modulation 
with  mf   =   0.6  where  the  f    =   1/10  Hz  and  fm   =   l/60Hz. 
As  can  be  seen,  the  two  spectra  are  almost  identical.   Both 
include  side  bands  centered  at  the  same  frequencies  about 
the  fundamental. 
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What  can  be  learned  from  this  discussion  that  may  be 
applicable  to  real  ocean  waves  is  that  side  bands  should  be 
expected  in  narrow  band  wave  records.   The  presence  of  side 
bands  is  probably  due  to  some  kind  of  amplitude  modulation 
rather  than  frequency  modulation  because  the  latter  inter- 
pretation implies  only  a  very  small  change  in  the  wave  periods 
(see  Figure  72)  and  a  constant  amplitude,  factors  normally  not 
expected  in  ocean  wave  records.   On  the  other  hand,  larger 
variability  in  the  periods  can  be  interpreted  as  the  conse- 
quence of  a  larger  m^,  and  necessarily  a  larger  number  of 
components  should  be  present  in  the  spectrum. 

Thus  it  seems  more  interesting  to  confirm  whether  the 
side  bands  that  can  be  identified  in  a  wave  spectrum  may  be 
associated  with  a  more  or  less  well-defined  amplitude  modu- 
lation which  appears  in  the  wave  record  as  wave-groups. 

B.   SIWEH  AND  WAVE  GROUPS 

When  our  research  on  the  side  bands  of  the  spectrum  had 
reached  this  point  a  very  interesting  report  by  Funke  and 
Mansard  (1979)  on  the  subject  of  synthesizing  wave  records 
was  brought  to  our  attention. 

For  the  purpose  of  isolating  wave  groups  in  a  wave  record, 
these  authors  note  that  "a  more  meaningful  description  of 
group  activity  could  be  obtained  by  computing  the  distribution 
of  wave  energy  along  the  time  axis  where  wave  energy  would 
be  defined  as  the  square  of  the  water  surface  elevation 
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averaged  over  a  period  which  is  a  function  of  the  peak 
frequency."   In  a  sense  this  operation  represents  a  process 
of  smoothing  the  time  series;  that  is,  if  we  look  in  the 
frequency  domain,  what  they  have  done  is  to  perform  a  filter- 
ing process  with  a  low-pass  band  filter  having  a  cut-off 
frequency  equal  to  the  peak  spectral  frequency. 

The  name  that  they  give  to  the  smoothed  wave  record,  the 
Smooth  Instantaneous  Wave  Energy  History  (SIWEH) ,  seems  to 
us  to  be  well  chosen  and  we  will  use  it  here. 

1.   SIWEH  Analysis 

The  Smooth  Instantaneous  Wave  Energy  History  (SIWEH) 
is  produced  by  operating  on  the  wave  record  using  the  Bartlett 
window,  which  is  a  triangular  weighting  function  defined  in 
mathematical  terms  by  the  following: 

SIWEH   =   E(t) 
T 


where:   E(t)   =   -==-  /  *  \  z  ( t  +  x)  .Q±  (t  )  di 

P 

*P 

for    T   <  t  -  T   -  T 
p        n    p 

and   Q, (t)   =   1  -  |t|/T   for  -T  <t<T 

J.  P  P      ir 

=   0  everywhere  else 
where:   Tn   =   duration  of  the  finite  wave  record 

T    =   peak  spectral  period. 
For  the  beginning  and  end  segments  of  the  record  we  have 
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~(T %    t)  I  P      Y\2(t   +    I)  .Q1(i)dx 

P  -^T     =     -t 

for        0    -    t   ^   Tp 

yT    -t 
11  +\2(t    +     T)    .Q1(T)dT 

e  T    =    -T 

p 

T      -    T       -    t    -    T 

The  product  of  the  SIWEH  operation  on  the  wave  record 
is  a  time  series  that  reflects  the  envelope  of  the  record, 
and  thereby  identifies  the  larger  waves  that  constitute  wave 
groups.   Examples  of  SIWEH  are  shown  in  Figures  4,  22,  23, 
33,  43,  53,  and  63  (upper  graphs) .   Another  comment  by  these 
authors  is  necessary  at  this  point.   In  their  analysis,  the 
authors  make  the  important  comment  that  their  work  was  done 
with  "paper  waves"  which  means  simulated  data.   Therefore  it 
was  considered  worthwhile  here  to  analyze  our  thoughts  about 
low  frequency  amplitude  modulation  utilizing  experiments  with 
real  wave  data. 

Also,  some  parameters  were  introduced  by  Funke  and 
Mansard  to  describe  wave  groups  that  had  not  been  evaluated 
using  real  wave  data.   We  took  the  opportunity  here  to  produce 
such  data  for  our  19  wave  records.   These  parameters  and  the 
results  obtained  are  described  in  the  following  sections. 
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2 .   SIWEH  Wave  Group  Parameters 

The  identification  of  the  low-frequency  components  of 
our  wave  records  was  done  by  the  application  of  the  already 
defined  SIWEH. 

This  analysis  of  wave  groups  is  based  on  the  concept 
expressed  by  Funke  and  Mansard  (19  79)  that  it  is  important  to 
characterize  at  least  four  parameters,  each  derived  from  the 
SIWEH,  in  order  to  fully  describe  wave  group  activity. 

The  first  parameter,  the  groupiness  factor,  indicates 
whether  there  are  wave  groups  present  and  in  what  degree. 
This  factor  is  defined  by 


"I —     n 


where:   E(t)   =   SIWEH 

In  words,  the  groupiness  factor  is  the  standard  deviation  of 

the  SIWEH  about  its  mean  and  normalized  with  respect  to  this 

mean. 

The  other  three  parameters  describe  the  average  group 
duration,  the  group  repetition  period,  and  the  variability  of 
the  group  repetition  period. 

To  calculate  the  group  duration  the  autocorrelation 
function  was  applied  to  the  SIWEH,  and  compared  with  the  group 
length  measured  in  the  following  manner: 
-  The  mean  energy  E  of  SIWEH  was  computed. 
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-  Using  E  as  a  threshold  the  intersections  between  SIWEH 
and  E  were  determined. 

-  The  slope  of  SIWEH  at  the  points  of  intersection  on  both 
sides  of  each  wave  group  was  projected  to  the  zero-energy 
line. 

-  These  intersections  with  the  zero-energy  line  define  the 
extent  of  the  group. 

For  the  group  repetition  period  and  its  variability 
the  SIWEH  spectral  density  function  was  computed  using  again 
four  "windows" ,  which  means  a  maximum  detectable  period  of 
256  seconds  with  a  resolution  of  1/256  Hz  and  eight  degrees 
of  freedom. 

The  results  of  these  computations  for  the  19  wave 
records  tabulated  by  wave  type  are  presented  in  Table  IV  and 
are  discussed  in  the  following  sections. 
3 .   Wave  Group  Properties 

Based  on  the  criterion  of  significant  wave  steepness 
the  19  wave  records  analyzed  were  divided  into  groups  by 
wave  type:   sea,  young  swell,  moderate  swell,  and  old  swell, 
according  to  the  definitions  given  earlier. 

For  illustrative  purposes,  one  wave  record  of  each 

type  was  chosen  as  follows: 

File  Number  579  Sea 

File  Number  600  Young  swell 

File  Number  446  Moderate  swell 

File  Number  677  Old  swell 

Other  illustrations  from  File  Numbers  228,  178,  and 
451  are  also  presented  as  examples  of  particular  character- 
istics that  will  be  pointed  out  later. 
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a.  Groupiness  Factor 

In  Table  IV  we  present  the  values  of  the  groupi- 
ness factor  found.   The  values  vary  widely,  but  show  a  weak 
relationship  to  the  wave  type  (see  Table  I)  with  sea  having 
the  highest  median  value  and  old  swell  the  lowest  value. 

Figure  4  (File  Number  228)  and  Figure  22  (File 
Number  451)  represent  wave  records  having  high  groupiness 
factors.   It  can  be  seen  that  the  higher  values  correspond 
to  wave  records  where  a  high  concentration  of  energy  appears 
after  relatively  long  periods  of  lower  activity  close  to  the 
mean  value.   It  is  interesting  to  note  that  we  did  not  find 
the  high  values  that  are  associated  with  the  synthesized  waves 
described  by  Funke  and  Mansard  (1979),  for  example  GF   =   1.35 

b.  Duration  of  Wave  Groups 

The  duration  of  wave  groups,  as  expected,  acquires 
some  meaning  only  with  swell-like  records.   Here  wave  groups 
are  usually  easily  identified  visually,  but  examination  of 
the  SIWEH  shows  that  a  typical  mean  duration  of  wave  groups 
is  not  easy  to  define.   In  sea-like  wave  records  the  wave 
groups  are  not  defined  and  the  interpretation  of  the  SIWEH 
has  no  apparent  meaning  other  than  to  signify  an  essentially 
random  variation  of  wave  heights. 

The  general  conclusion  from  these  19  wave  records 
analyzed  is  that  the  zero-crossing  time  lag  from  the  auto- 
correlation graphs  (see  Figures  16,  30,  40,  50,  and  60)  does 


48 


not  describe  the  mean  duration  of  wave  groups,  but  that  a 
rough  estimative  for  this  parameter  can  be  obtained  from  the 
SIWEH . 

c.   Group  Repetition  Period  and  Variability 

The  spectrum  of  SIWEH  for  swell  records  gives  an 
indication  of  the  periodicity  of  wave  groups  (Figures  51  and 
61) .   The  values  found  and  presented  in  Table  IV  are  condi- 
tioned by  our  256  second  "window"  and  the  consequent  resolu- 
tion of  1/256  Hz.   This  implies  that  the  peak  frequency  should 
be  a  multiple  of  the  fundamental  frequency  1/256  Hz  such  as 
that  shown  in  Figure  20,  and  this  explains  the  periods 
presented.   It  may  be  noted  that  these  periods  generally 
coincide  with  the  periods  of  surf  beat  commonly  seen  in 
coastal  wave  records. 

As  discussed  above,  the  presence  of  wave  groups 
may  be  indicated  by  the  presence  of  side  bands  that  appear  in 
the  power  density  spectrum  according  to  the  theory  of  ampli- 
tude modulation.   Thus  Figure  59  of  File  Number  6  77  shows  two 
side  bands  present  in  the  PDS  that  are  symmetrically  distri- 
buted about  the  peak  spectral  period  of  15  seconds.   This 
low-frequency  signal,  also  visualized  in  the  envelope  of  the 
wave  record,  should  have  a  spectral  profile  identical  to  the 
side  bands.   The  PDS  of  the  SIWEH,  Figure  61,  confirms  this 
characteristic.   A  similar  situation  is  shown  in  Figures  10 
and  17  of  File  Number  228. 
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We  recognize  that  it  is  not  evident  from  the  PDS 
of  the  sea-surface  elevation  which  spectral  component  of  the 
side  bands  represents  the  peak  frequency  such  as  was  produced 
with  the  simulated  waves  (Figures  68  and  69) .   Thus  it  is 
important  in  defining  the  group  repetition  period  to  identify 
this  low  frequency  signal,  and  the  application  of  the  SIWEH 
is  one  process  for  doing  this. 

Thus,  as  the  PDS  of  the  SIWEH  of  the  simulated 
waves  (Figure  70)  shows  the  group  repetition  period  of  the 
simulated  wave  groups,  the  PDS  of  the  SIWEH  of  File  Number 
6  77  (Figure  61)  presents  a  group  repetition  period  that  agrees 
with  the  occurrence  of  wave  groups  visualized  on  Figure  53 
(upper  graph) . 

The  variability  of  the  wave  group  period  can  be 
seen  in  the  bandwidth  of  the  spectral  representation  of  the 
SIWEH.   A  broad  bandwidth  is  an  indication  of  large  varia- 
bility and  a  narrow  one  an  indication  of  small  variability. 
In  Table  IV  we  categorize  the  variability  only  as  large  or 
small  because  of  the  difficulty  of  establishing  some  quanti- 
fied parameter.   It  may  be  seen  that  the  variability  of  wave 
group  period  tends  to  be  large  with  young  swell  and  small  with 
old  swell.   The  form  of  the  PDS  of  the  SIWEH  presented  in 
Figures  17,  31,  41,  51,  and  61  is  instructive  when  compared 
with  the  simulated  example  presented  in  Figure  70  where 
theoretically  the  variability  is  zero. 
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V.   CONCLUSIONS 

This  thesis  presents  the  results  of  experiments  in  the 
analysis  of  ocean  wave  records.   For  this  purpose  19  digitized 
records  obtained  from  a  wave  buoy  in  60  meters  depth  off  the 
California  coast  were  analyzed.   Two  general  approaches  were 
taken:   (1)  examination  of  the  interrelationship  of  the 
properties,  H  and  T,  of  the  individual  waves  in  a  wave  record, 
and  (2)  study  of  the  properties  of  the  larger  sets  of  waves 
in  a  wave  record,  termed  wave  groups,  and  consideration  of 
wave  groups  as  a  manifestation  of  amplitude  modulation.   The 
principal  findings  are  as  follows. 

Wave  heights  and  periods,  determined  using  the  zero-upcross 
method  and  treated  as  a  function  of  H   versus  1/T  can  be 
manipulated  to  produce  a  cumulative  pseudo-energy  distribution 
that  approximates  the  cumulative  true  power  density.   Thus  it 
is  possible  to  reasonably  approximate  some  important  proper- 
ties of  the  PDS  by  the  use  of  the  H  and  T  values  obtained 
from  a  wave  record. 

Sequences  of  H  and  T  for  each  wave  record  were  produced 
graphically  in  order  to  provide  different  ways  of  looking  at 
the  wave  records  visually.   The  presence  of  groups  of 
successive  large  waves  is  readily  seen  in  the  wave  height 
sequence  plots.   Examination  of  period  sequences  plots  show 
these  larger  waves  to  have  almost  constant  and  longer  periods. 


51 


Plots  of  H  versus  T  were  also  produced.   Those  give  an 
indication  of  how  the  individual  wave  steepness  is  scattered 
about  the  significant  record  steepness. 

As  a  manifestation  of  amplitude  modulation,  the  side 
bands  in  the  PDS  of  the  sea  surface  elevation  give  for  narrow 
band  spectra  (swell  type  waves)  an  indication  of  wave  group 
activity  and  periodicity. 

The  Smooth  Instantaneous  Wave  Energy  History  (SIWEH) , 
developed  by  Funke  and  Mansard  (19  79)  from  work  with 
synthesized  waves,  and  the  power  density  spectrum  of  the 
SIWEH  were  found  to  provide  a  good  additional  means  of 
characterizing  wave  groups  in  narrow  band  ocean  wave  records 
and  to  give  information  on  the  frequency  of  occurrence, 
variability,  and  duration  of  wave  groups.   They  proved  less 
useful  for  studying  wide  band  spectra  (sea  waves) . 

The  groupiness  factor  derived  from  the  SIWEH  was  found 
to  be  a  fair  descriptor  of  the  degree  of  development  of  wave 
groups.   The  application  of  the  autocorrelation  function  to 
the  SIWEH  was  found  to  give  inconclusive  results  about  the 
duration  of  wave  groups. 
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TABLE  II 

TOTAL  ENERGY  PROPORTIONALITY  FACTOR 
BETWEEN  PSEUDO  AND  TRUE  SPECTRA 


N 


File     Number 
Number  Waves-N 


2 rr  2  i „J1\         Statistical 

.     :  H;    cm    tt   •      i       2- 
i=l   -1         Variance  (cm  . 


Factor  of 
Proportionality 


579 

151 

5238301 

5102 

6  .80 

683 

144 

6639113 

6384 

7.22 

336 

19  2 

2091401 

1768 

6.16 

686 

156 

3507162 

3  29  4 

6.83 

547 

169 

4813661 

4242 

6.71 

228 

123 

2863876 

3242 

7.18 

339 

144 

3129877 

3182 

6.83 

023 

146 

1659909 

1656 

6  .87 

600 

10  8 

5161522 

6949 

6.88 

611 

118 

5067237 

616  8 

6.96 

536 

119 

3847429 

4749 

6  .81 

612 

116 

4124692 

5333 

6  .67 

451 

128 

3369897 

4307 

6  .11 

410 

146 

4494133 

4597 

6.70 

607 

128 

463045 

528 

6.85 

446 

99 

2178222 

3  29  2 

6.68 

370 

192 

2042953 

1791 

5.94 

677 

97 

925187 

1360 

7.01 

178 

138 

1361507 

1632 

6.05 

Mean    6.70 

Standard  Deviation 

0  .37 
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FIG  2.  EXAMPLE  OF  CEDN  TABULAR  SPECTRA  PRESENTATION. 
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